Preprint typeset in JHEP style - PAPER VERSION 



Cavendish-HEP-05/03 



On the quark distribution in an on-shell heavy quark 
and its all-order relations with the perturbative 
fragmentation function 



Einan Gardi 

Cavendish Laboratory, University of Cambridge 
Madingley Road, Cambridge, CBS OHE, UK 



Abstract: I present new results on the quark distribution in an on-shell heavy quark 
in perturbative QCD and explore its all-order relations with heavy-quark fragmentation. 
I first compute the momentum distribution function to all orders in the large-/3o limit and 
show that it is identical to the perturbative heavy-quark fragmentation function in the same 
approximation. I then analyze the Sudakov limit of the distribution and the fragmentation 
functions using Wilson lines and prove that the corresponding Sudakov exponents in the 
non-Abelian theory are the same to any logarithmic accuracy. The anomalous dimension is 
then determined to two-loop order, corresponding to next-to-next-to-leading logarithmic 
accuracy in the exponent, in two ways: the first by extracting the singular terms from a 
recent calculation of the fragmentation function and the second by performing the two-loop 
Wilson-line calculation in configuration space. I find perfect agreement between the two. 
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1. Introduction 

Relations between space-Uke processes that probe the parton distribution function and 
time-fike processes involving fragmentation have been known since long. This includes 
the Drell-Yan-Levy relation between deep-inelastic structure functions and single-particle 
inclusive cross section in e+e~ annihilation [1], as well as the Gribov-Lipatov reciprocity 
relation [2], relating the space-like and timc-likc splitting functions that determine the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [3] evolution of parton densities and 
fragmentation functions, respectively In dimensional regularization, the latter relation is 
violated beyond the leading order. 
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In this paper we find and explore new all-order pcrturbative relations between distri- 
bution and fragmentation functions, which are specific to heavy quarks^. When dealing 
with heavy quarks one can define and compute perturbative distribution and fragmentation 
functions, replacing the hadronic states by ones composed of on-shell quarks and gluons. 
Owing to collinear singularities this is not possible for light quarks, where only the ultravi- 
olet (DGLAP) evolution of these functions can be computed in perturbation theory, while 
the actual matrix elements, which set the initial condition for the evolution, can only be 
defined non-perturbatively. For heavy quarks the mass effectively cuts off collinear radia- 
tion and, if m S> A, it provides a hard scale for the coupling. This makes the perturbative 
initial condition for DGLAP evolution of the distribution and fragmentation functions well 
defined to any order in perturbation theory. These pcrturbative initial conditions are the 
subject of this paper. Of course, these functions still differ from their non-perturbative 
counterparts by power corrections, which wc shall briefiy discuss as well. 

The non-perturbative process-independent definitions of distribution and fragmenta- 
tion functions [4,5] — see Eqs. (2.1) and (3.30), respectively — apply to both light and 
heavy quarks. These field-theoretic definitions are based on the Fourier transform of certain 
non-local matrix elements on the lightcone. The distribution function f{x; jj) measures the 
longitudinal momentum distribution of a qiiark in a hadron, where the quark momentum 
fraction is x G [0, 1]. This function can be interpreted as a probability distribution, barring 
the fact that it requires renormalization. Similarly, the fragmentation function d{x\ ji) mea- 
sures the probability distribution of producing a hadron with fraction x of the momentum 
originally carried by the quark. Also here x G [0, 1]. Using the heavy quark expansion [6] 
one can show that for /x ~ m both /(x; /x) and d{x; /x) peak at 1 — a; = 0{A/m). 

The momentum distribution function of a heavy quark in a heavy meson plays a central 
role in the calculation of inclusive decay spectra such as B — > Xg^ and semileptonic B 
decays [7 14] . This distribution essentially determines the shape of inclusive decay spectra 
in the experimentally-important endpoint region. In this kinematic domain the invari- 
ant mass of the produced hadronic system is small compared to the quark mass and the 
0{A/M) fraction of the momentum which is carried by the light degrees of freedom in the 
meson becomes important. In this context it is convenient to analyze the quark distribution 
function in the infinite-mass limit, where it is often called the "Shape Function". 

Although the quark distribution in a heavy meson is a non-perturbative distribution, it 
has an important perturbative ingredient. Being completely inclusive, this distribution can 
be approximated by the quark distribution in an on-shell heavy quark fprix; //), while non- 
perturbative effects are treated as corrections. This idea is theoretically appealing, however, 
it is not easy to achieve. First, the pcrturbative description of the quark distribution in an 
on-shell heavy quark can only be reliable if Sudakov logarithms are resummed, resummation 
that takes the form of exponentiation in moment space. Moreover, running-coupling effects 
are significant and lead to infrared renormalons in these moments. The corresponding 
ambiguities are only resolved at the non-perturbative level, making the resummation of 
running-coupling effects necessary for systematic separation between perturbative and non- 

^We do not consider in this paper deep inelastic structure functions. 
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perturbative contributions. The Dressed Gluon Exponentiation (DGE) approach [14-17,30] 
incorporates renormalon resummation in the calculation of the Sudakov exponent. This 
opens up the way for consistently using the resummed quark distribution in an on-shell 
quark as the baseline for describing the quark distribution in a heavy meson. 

In a completely different application of QCD, namely heavy-quark production in hard 
scattering processes, one encounters the heavy-quark fragmentation function, d{x; n). Also 
here the resummed perturbative initial condition for the evolution, dpT{x;n), provides a 
baseline for a systematic description of the non-perturbative object [27-31]. DGE-based 
predictions [30], which involve just one or two non-perturbative power corrections, have 
proven successful in describing the inclusive b production cross section measured at LEP. 

In this paper wc compute the quark distribution in an on-shell heavy quark, fprix; /x), 
in two different limits. In each case a certain gauge— invariant set of radiative corrections 
is controlled to all orders. The two limits are: 

• The single-dressed-gluon approximation (Sec. 2). Taking the large-/3o limit — or, 
formally, the large-A^j limit — we perform an all-order resummation of running cou- 
pling effects (renormalons) . This resummation has two applications [18-20]: first, it 
improves the leading-order calculation incorporating BLM-type radiative corrections 
to all orders [21,22], and second, it probes the structure of power corrections. 

• The X — > 1 Sudakov limit (Sec. 3). Taking the large-x limit we keep only singular 
terms in fprix; //), which build up the Sudakov exponent. Here the calculation is done 
to two-loop order, corresponding to next-to-next-to-leading logarithmic (NNLL) 
accuracy in the exponent. More generally, following the work of Korchcmsky and 
Marchesini in Ref . [23] , we show that the singular x — > 1 terms are fully captured by 
the distribution defined in the m — > oo limit, the so-called "Shape Function". This 
means that the exponent is computable, to all orders, through the renormalization 
of a Il-shaped Wilson-line operator^ with two antiparallel timelike rays, connected 
by a lightlike segment. 

Interestingly, in both limits we find that exactly the same resummed expressions hold for 

the perturbative heavy-quark fragmentation function [29-31]. In general though, the two 
functions do differ; for example, their DGLAP evolution away from the large-x limit starts 
differing at two-loops [2]. 

2. Quairk distribution in an on-shell heavy quEirk computed with a single 
dressed gluon 

2.1 Definition and calculation of the gluon emission cut 

We define the parton distribution [4] in a heavy meson \H{p)) in the standard way, by 

f{x;n) = r '^e-'-P^y- (F(p)|*(y)$,(0,2/)7+*(0)|iJ(p))^, (2.1) 

^Recall that not all Sudakov anomalous dimensions can be computed in the Eikonal approximation. For 

example, the jet function controlling the large- a; limit of deep-inelastic structure functions [16,24-26], which 
is sensitive to coUinear radiation from a light quark, cannot be reproduced in this approximation. 
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where ^ is introduced as an ultraviolet renormalization scale for the operator, y is a light- 
like vector in the "— " direction, ^y(Q,y) is a path-ordered exponential in this direction 
connecting the points y and 0, i.e. 



^y(0> y) = 'P exp (^ig j dz^j^A^ 



(^) 



Some properties of $j^(0, y) are collected in Appendix A. 

In order to compute the matrix element it is convenient use the lightcone gauge, where 

E^'^^^'^^^-g/^^ + ^'V^'^^" - (2-2) 

Then, at one loop there is only one diagram: the box. The calculation, based on of Eq. (2.1) 
is summarized in Appendix A. In 4 — 2e dimensions the result for the not-yet-renormalized 
quark distribution is: 

#,-'^)(x;/.) = (5(l-x)[l + 0(«.)]+(^^J duT{u)^-^j S(x,«,e) (2.3) 



with 



B{x, u, e) = I|±^ ef " (1 - x)-^-^ V 



(2.4) 



where we used the Borel-modified gluon propagator l/(— A;^) — ^ l/(— fc^)^^" in order to 
resum running-coupling effects of all orders. With one-loop running coupling T(n) = 1. To 
go beyond one-loop running we use the scheme-invariant Borel representation [32] where 
T(u) is the Laplace transform of the 't Hooft coupling: 



with 5 = Pi/Pl where po = %Ca - ^Nf and pi = - gCpiV/ - ^C^iV/. We define 

A in MS. 

Out of the terms in the square brackets in Eq. (2.4) only (with a coefficient 1) 
comes from the Axial gauge part of the propagator in Eq. (2.2), which is proportional to 
Xjk ■ y. Through the y~ and the momentum integration this singularity of the propagator 
is converted into a 1/(1 — a;) singularity, k ■ y/p ■ y — > (1 — x). The square brackets are 
a generalization of the splitting function. It is noted that Eq. (2.3) is identical to the 
perturbative fragmentation function [29] computed in the same approximation in Ref. [30] 
using a different technique. 

2.2 Moments and evolution 

In Eq. (2.3) we computed only the real gluon emission contribution. Clearly this result is 
singular at x — > 1. Virtual corrections generate 0{as) divergent contributions, which are 
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proportional to (5(1 — x). Owing to the inclusive nature of Eq. (2.1) infrared singularities 
cancel out. However, there are also ultraviolet singularities. We therefore proceed by 
first taking a logarithmic derivative with respect to the mass to remove the ultraviolet 
divergence. Then we can go to four dimensions setting e = 0. Next, we go to moment 
space, 

dxx^-^f{x-iJL), (2.6) 



Fn 



where wc can use the fact that the first moment Fn=i corresponds to the total number of 
quarks in the quark, a conserved current, so it is -pAr=i = 1 for any mass and the derivative 
must vanish identically. This allows us to reconstruct the virtual (iV-independent) terms 
missing in Eq. (2.3) out of the real-emission (AT-dependent) ones. We get: 



dlnFN 
din rn? 



Cf 



r /"OO 

1^ / duT{u) 
'0 Jo 



A2 
m 



2 , Bf{N,u), 



(2.7) 



with 
B 



■f{N, u)= j dx [x^-^ - 1) B{x, M, e = 0) = 
Jo 

T{N + l + u) r(2 + u) 



(1 - u)r(-2^x) 



r(Ar+l-n) r(2-n) 



1 



+-{l + u)T{2-2u 



J r(A- + ») r(i + tx) y 

' \ V{N + 2-u) r(3-n)y_ 



(2.8) 

+ O{u/f3o). 



Next, we would like to integrate over Inm^ to recover F^. However, this brings back 
the ultraviolet divergence which takes the form of a l/u singularity. Therefore, we must 
perform ultraviolet subtraction. The result is: 



TTiPT 



Cf du , , 
Po Jo u 
+ Oil/Pi), 



4)" 



Bf{N,u) + 



m 



2\ " 



Sap(A^,w) 



(2.9) 



where ^Ui? is a factorization scale and Bap{N,u) is the (non-singlet) Altarelli-Parisi evolu- 
tion kernel. 



Bap{N,u) = ^ 



n=0 



(2.10) 



The leading order coefficient in (2.10) is renormalization-scheme invariant, and it equals 
to the u = limit of (2.8), ensuring the cancellation of the 1/u singularity in (2.9): 



7o(A^) = Si(7V)-^ + l 



1 



N + 1 



1 

N 



where Sk{N) = J2f=i so Si{N) = '^{N+l)+'yE- In the MS factorization scheme 7„(iV) 
are known to NNLO (n = 2) in full [33] and to all orders in the large-/3o limit [34], see 
Appendix C. 
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We find that the final result for the moments of the quark distribution in the large-/3o 
limit, given by Eq. (2.9), is identical to the one corresponding to the heavy-quark fragmen- 
tation in the same approximation, see Eqs. (38) and (39) in Ref. [30]. The expansion of 
FJ^ijn; fip) in powers of Og can be readily obtained from Eq. (2.9) by expanding the terms 
in the square brackets in powers of u and integrating order by order. The result, to 0(0^), 
appears in Eq. (43) of Ref. [30]. 

The infrarcd-renormalon structure can be read off Eq. (2.8): there are simple poles at 
all integer and half integer values of u, except for u = 1. Note that all these singularities 
are absent in the real-emission expression in x space, where one finds just an upper limit 
of the form {1 — x) A/m on the values of x for which the Borel integral exists. These 
renormalons are all associated with the integration over x for x — > 1, and within the 
scheme-invariant formulation of the Borel transform (where /?i -terms are factored out) they 
are expected to remain simple poles in the full theory. As already noted in Refs. [14,30], 
the residues in Eq. (2.8) depend on A'^ such that at large A'^ ambiguities appear as integer 
powers of (NA/m). These parametrically-enhanced power terms are related to renormalon 
ambiguities in the Sudakov exponent, and they will be discussed further in the next section. 

3. Sudeikov logs in the queirk distribution function 
3.1 Definitions and results 

In the large- AT limit Sudakov logarithms exponentiate as follows: 

F^^(m; fip) = H{m; hf) 5iv(m; hf) + 0{1/N), (3.1) 



where 



n=l 
1 



TT 

_ 1 



'71+ 1 



k=l 



(3.2) 



exp < I dx 

/o l-a? 



2dfj, 



(l—x)m A* 



Here^ A and V are Sudakov anomalous dimensions having the following perturbative ex- 
pansions in the MS renormalization scheme: 

AM.)) - 1 ^ (^)" - ^ E"? (^)" . 

n=l ^ ^ n=l ^ ^ 

= f Z,-(<M)" ^ f f;.-(*i&))", (3.3, 

n=l ^ / yo \ J 



^Notations in the literature vary; for example in Ref. [29], Eq. (69), our A and V are denoted by A 
and —H, respectively. 
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where the normahzation^ is such that ai = di = 1. The coefficients 02*^ and af^ are 
known based on calculations in deep inelastic scattering (see e.g. [33]). They are given in 
Eq. (C.3). df^ is given by: 

^-l^^dC:.-^^-'^)- (3.4) 



9 Po 12 18 

We shall explain how it was determined in Sec. 3.5. 

A{as{n)) is the universal cusp anomalous dimension which generates the double logs. 
It is well known that this object can be defined through the renormalization of a Wilson 
line with a cusp [23,35-40]. P(as(/u)) is another anomalous dimension, which describes 
soft radiation from the heavy quark and it is specific to the heavy-quark distribution 
function. In Sec. 3.4 we shall prove that the same anomalous dimension controls the large- 
N limit of the heavy-quark fragmentation function. P(as(//)) too can be defined and 
computed using a Wilson-line operator. In contrast with the cusp anomalous dimension, 
it is associated with a specific Il-shaped configuration with a finite light-like section to 
which two antiparallcl infinite time-like rays attach. This object was first analyzed by 
Korchemsky and Marchcsini in Rcf. [23]. 

The single-dressed-gluon calculation of the previous section has just leading logarith- 
mic accuracy: Eq. (2.8) (or Eq. (3.10) below) fixes Cn,k=n+i for any n, but gives only the 
part of the Cn^k<n which is leading in (3q. In the full theory there are additional terms 
in the exponent which have different color factors. While A (the cusp anomalous dimen- 
sion) contains both non-Abelian Ca/Pq and Abelian Cp/Po terms, T) has only non-Abclian 
ones. These coefficients can be determined by computing the anomalous dimensions order 
by order. Next -to-leading logarithmic accuracy (Cn,fc>n) requires fixing 02 and di. The 
state of the art is next-to-next-to-leading logarithmic accuracy {Cn,k>n-i) which requires 
knowing also 03 and ^2 . Note that the calculation of the exponent to this formal accuracy 
should be complemented by full next-to-next-to-leading order calculations of the hard 
coefficient function. These are yet unavailable for radiative or semileptonic decay spectra. 

Switching to the scheme invariant Borel representation of the two anomalous dimen- 
sions, 

A{asii2)) = -f duTiu) ( — ) BAiu), 
Po Jo / 

V{as{ix)) = -f dunu)( — ] Bj,{u), (3.5) 
Po Jo J 

the moments of the distribution function become: 

F^^(m; hf) = H{m; iip) 5iv(m; lip) + C»(l/iV), (3.6) 



''We shall mostly use the notation a„ and d„ for the coefficients. This is convenient for comparison with 
the large i9o limit and with the Borel formulation but it does not imply any additional approximation: On 

and dn contain all color factors. 
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with^ 



f t>3 , , \\n n+l >. 
n=l V ^ / fe=i J 



(3.7) 



expi^/ -r(u)(-^ 



S5(n)r(-2u) (Ar2" - 1) + ( 11^ ) 5^(u)lniV 



m 



Po Jo u 

where Cn,k=i depend on ln/Xi?/m while Cn,k>i are numbers, and where 

Bsiu) = BAiu)-uBv{u). 



(3.8) 



These expressions are completely general. 

Based on Eq. (2.8) we find that the hard function that incorporates the finite terms 
at N — ^ oo is given by: 



H{m\iiF) = 1 + 



CfO-s 



TT 



3 \ ^ TT^ 2 

-T + 7B In — + 1 - -7- + 7E - 7e 

4 / 6 



+ 0{a^), (3.9) 



and that in the large-/3o limit the Borel transform of the anomalous dimension Bs{u) is 
given by 



Bs{u) = e3«(l-u) + C»(l//3o). 



(3.10) 



The information contained in Eqs. (3.10), (3.8) and (C.4) can readily be translated into 
values of the coefficients in Eq. (3.3) in the large-/3o limit. One gets: 



_ 1 

large 0a Tl 



large /3o 



3n 



(3.11) 



The first few orders are summarized in Tabic 1. Note that the convergence of the cusp 
anomalous dimension A is much faster® than that of T>. This property presumably persists 
in the full theory. 

Expanding the Borel integral and the exponential in Eq. (3.7) one obtains the log- 
enhanced terms, order by order in the coupling. In contrast with Eq. (3.2), Eq. (3.7) 
incorporates renormalon resummation in the Sudakov exponent (DGE). It exposes the fact 
that the exponent contains infrared renormalon ambiguities: poles in T{—2u) along the 
integration path, which generate power-like ambiguities ~ {NA/my , where j are integers. 
These parametrically-large ambiguities are indicative of corresponding non-perturbative 
power terms. As shown in Ref. [14] (see also [39]), the u = ^ renormalon is related to the 
leading infrared renormalon in the pole mass [41,42], while higher powers are associated 



^Although both Eq. (3.2) and Eq. (3.7) axe normalized such that the N = 1 moment is identically 

unity, the integration over x in Eq. (3.2) generates, in addition to the relevant InA'^ terms, some finite 
terms along with terms that vanish as powers of For this reason the functions in Eq. (3.1) have a 

tilde distinguishing them from those of Eq. (3.6). As far as the logarithms are concerned the exponent in 
Eq. (3.7) is identical to that of Eq. (3.2). 

®In fact it is probably faster than all other Sudakov anomalous dimension. See e.g. table 1 in [26]. 
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n 




larRe 00 




large 0o 


ratio 


i 


1 




1 




i. 


2 
3 
4 
5 
6 


5 
3 

_ 1 

i-2C3 

30 " 3 3 

^-6C5 + ^7r^ + iC3 


1 
9 

109 2 a 

81 3 ''3 
1 „4 I 212 5 A 

120 " ""^ 81 6 ''S 

6331 6 a I 1 _4 1 2 a 

1215 5''5 -1- 90 " -r- i^i^3 

5 a 1 _i_ 1 ^6 I 20191 _L 1 a2 1a 
3^5 540 " "T 567 2187 '^3^3 9^^ 


15. 

-0.6126 
-0.8531 

-0.2099 
0.0347 



Table 1: The large-/3o part of the coefficients ajf® and djf® in Eq. (3.3) and their ratios. 



with other local matrix elements which constitute the non-perturbative quark distribution 
function [8] or "Shape Function" [7]. 

In the DGE approach the Borel integration is performed using the Cauchy Principal 
Value prescription, making an explicit power-like separation between perturbative and non- 
perturbative contributions to the exponent. Any information on the anomalous dimensions 
from fixed-order calculations can be used in Eq. (3.7). Based on Eqs. (C.3) and (3.4), we 
have: 

BA{n) = l+(l + C2]'^+(-l+ c^) '4 + 0{u% (3.12) 



Bviu) = 1 + 



^3 J 11 \ 3 J 2! 
9 (3o V4 12 18 



^ + 0{u% (3.13) 



where c„ represent contributions that are subleading in Pq; C2 and C3 are given in Ap- 
pendix C. Note that the relations of c„ for n > 3 with ajf® involves the coefficients of the 
P function. The same is true for the 0{v?) and higher order coefficients in Bx>{u). Fi- 
nally, note that since the analytic dependence on u is not known beyond the large-/3o limit, 
there remains some uncertainty in evaluating the u-integral in Eq. (3.7) with Eqs. (3.12) 
and (3.13). This issue is addressed in detail in a forthcoming publication [17]. 

3.2 Strictly factorized form and Wilson lines 

The factorization-scale dependence of the Sudakov factor of Eq. (3.7) (or Eq. (3.2)) is 

governed by the cusp anomalous dimension alone while its dependence on the soft scale 
m/N is different. Furthermore, by convention this factor is normalized such that its first 
moment (A'^ = 1) is unity for any fip; this normalization is natural since the full quark 
distribution function indeed obeys F]\f=i = 1. These constraints can be satisfied only if 
the exponent acquires some dependence on the hard scale m. However, such dependence 
cannot be consistent with the Wilson-line operator definition [10,23,35,38], nor with the 
effective field theory approach [11-13,43,44]. In the latter definitions the Sudakov factor 
cannot depend on m, but only on the soft scale m/N and on the factorization scale, the 
renormalization scale of the operator. In order to convert Eq. (3.7) to a strictly factorized 
form we reshuffle radiative corrections that depend only on the hard scales into the hard 
function, writing 

Fj^^{m-iiF) = H{m;iiF)SN{m;iiF)+ 0{l/N)=n{m;iiF)S{NiiF/m)+ 0{l/N) (3.14) 
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where the first expression represents the definition of Sec. 3.1 whereas the second assumes 
strict factorization, as in Ref. [23]. Although we do not write it exphcitly, both 7i and S 
depends also on as{n), where the renormalization scale of the coupling, is not necessarily 
equal to the factorization scale fxp. In practice we shall be using the standard dimensional 
regularization with n defined in the MS scheme for both the operator and the coupling. 
We can now write the Borel representation of »S(iV/x/m) as follows: 

S{Nfi/m)=e^plf2(^y X]Cn,fcln'(iV/x/m)|= (3.15) 



Cf f°° du ^, , /A' 



exp / —T{u)[ — 

/3o Jo u V/^ 



2 \ " 



Bs{u)T{-2u) {{Nfi/my'' - l) + Ba{u) ln(iV/i/m) 



Bs{u) and Bj,{u) arc the same Borel functions appearing in Eq. (3.7); clearly Eq. (3.15) 
differs from Eq. (3.7) just by A'"-independent terms, and it is therefore consistent with 
Eq. (3.14). Note that in contrast with Eq. (3.7) the normalization (the N = 1 moment) 
of the strictly-factorized Sudakov factor of Eq. (3.15) strongly depends on Inix/m to any 
order in perturbation theory. 

Since there are no collinear singularities from the heavy-quark line, Sudakov logarithms 
in S{Nji/m) are related to soft gluons and to gluons that are collinear to the lightcone 
direction y in Eq. (2.1). One therefore expects that the log-enhanced terms of the quark 
distribution function Eq. (2.1) would not change if the heavy-quark lines are replaced by 
time-like Wilson lines. The Wilson line describes the Eikonal interaction of the heavy 
quark with soft gluons. This is equivalent to taking the quark mass to infinity. 

In order to convert Eq. (2.1) to a Wilson-line operator definition [23] one applies the 
Eikonal approximation replacing the interacting quark field '^{z) by a free heavy quark 
field il){z) multiplied by a Wilson line in the direction p which extends to infinity, i.e. 

^{z) $p(oo, z)il^{z) ■ *(z) ijj{z)^-p{z, oo) (3.16) 

with ^p{zi,Z2) defined as in Appendix A. Next, one uses the fact that the free fields 
annihilate the external quark states to convert the matrix element to one in the vacuum. 
Factoring out the Dirac structure one obtains: 

{h{p)\^iy)%iO,y)j+^iO)\hip))^-^ 2p+e^y-p^ (0| $_p(y, oo)a>,(0, y)$p(oo, 0) |0)^ . (3.17) 
The X — ^ 1 singular terms in the quark distribution function are therefore: 

/pt(x;m)L_i - fj^^"^''^'""^ (0|$_p(2/,oo)$,(0,y)$^(oo,0)|0)^, (3.18) 
Following Ref. [23] we define W[Cs\ by 

W[Cs]{ip ■ yii/m, as{^l)) = (0| $_p(2/, oo)$^(0, y)$p(oo, 0) |0)^ ; (3.19) 
it is shown in Fig. 1. W[Cs\{ip-yiJ'/m, as {/J,)) is an analytic function except on the negative 
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y 




Figure 1: Minkowski space-time picture (vertical axis as time and horizontal axis as x^) of the 
Wilson-line configuration W[Cs]{ip- yn/m) of Eq. (3.19) representing the quark distribution func- 
tion in an on-shell quark in the infinite mass limit (in the rest frame of this quark) . The two figures 
describe the situation when y~ is positive (l.h.s) or negative (r.h.s), where path-ordering on the 
lightlike segment I2 from to y corresponds to time-ordering and anti-time-ordering, respectively. 



real axis of ip ■ yji/m. We note that, since fpT{x; IJ,) in Eq. (3.18) is real, sign inversion of 
p ■ y is equivalent to complex conjugation, i.e. 

W[Cs]{ip ■ Vli/m, asifi)) = W[Csn-ip ■ Vli/m, «,(//)). (3.20) 

Up to 0{1/N) corrections the Fourier integral in Eq. (2.1) and the Mellin integral 
(Eq. (2.6)) amount to replacing ip~^y~ by N [23,39]. This can be shown by first extending 
the X integral to —00 relying on the fact that the contribution is only from the x near 1 
region, then performing the x integral getting —l/[ip^y~ — (N — 1)], and finally evaluating 
the y~ integral by residue, closing the contour through the lower half p- y plane. The result 
is 

F^^(m;M) = H^{m;i^) x W[Cs]{Ni^/m,as{fi)) + 0{1/N), (3.21) 

where Hy^{m; fi) is a hard coefficient function. It accounts for finite terms at N — > 00 
that are not captured by the Eikonal approximation of Eq. (3.17). The NLO expression 
for this function is given in Eq. (3.24) below. WfCg] is therefore related to S of Eq. (3.15) 

by 

where the additional normalization is needed to remove all A'^-independent pieces, which 
are absent in Eq. (3.15) by construction. 
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The Wilson line operator VF[C5] was analyzed in detail in Ref. [23], where its evo- 
lution equation was derived and a two-loop calculation of the corresponding anomalous 
dimensions was performed. Here we repeated this calculation finding: 



hlW[Cs]{N^i/m,as{|l)) = Cp 
+ Cf 



TT 



18 9 



+ 



55 9 ^ 77r2 

h - Cs 

108 4 18 



(3.23) 



TT 



+ 



where L = In Njl/m with p, = fie^^ where /i is defined in the MS scheme. 

The same logarithmic terms can be obtained from Eq. (3.15) by expressing it in terms 
of L, 



Cf 



Po Jo u 



2 \ " 



Bsiu)T{-2u) (e2(^-T^)" - l) +BAiu)iL-jE) 



substituting the anomalous dimension coefficients of Eqs. (3.8), (3.12) and (3.13) and ex- 
panding. The A/^-independent terms, such as — ^7r^CirQ;s(/i)/7r, are of course different, see 
Eq. (3.22). For H^{m;n) in Eq. (3.21) we find: 



H^{m-ii) = 1 + Cf- 



TT 



, 9 m 1 , m TT 
In In h — + 1 



IX 24 



+ •••• 



(3.24) 



Our result in Eq. (3.23) agrees with that of Ref. [23] except for the rational number 
in the coefficient of Ca in the single log term at order , where we find — ^ rather than 
~W appearing in Eq. (3.6) of Ref. [23]. This discrepancy is due to an error in Ref. [23] 
in the evaluation of one of the diagrams — see Sec. 3.5 and Appendix D below. 



3.3 Evolution equation 

From Eq. (3.15) above one can be obtain an evolution equation by taking a full logarithmic 
derivative with respect to the scale. We obtain: 



dlnS{Nn/m) ^ d\nS{Nn/m) das{n) 



d In fi 



dag (/x) dln/j. 
d\nS{Nfi/m) 
In jji 



(3.25) 



2A{as{lx)) In— ^ + Vs{as{ix)) 



m 



where, as above, jl = fx e^^ , and the anomalous dimension F^ can be related order by order 
to the two anomalous dimensions V and A defined above: 



Cf f°° 
TsMi^)) = -3- / duT{i 
Po Jo 



2 \ f 



F(l - 2u)Bviu) + 



1 + 2n7B - r(l - 2u) 



u 



Ba{u) 
(3.26) 
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For its expansion we obtain: 



TT 



+ /3o 



Cf 



TT 



TT 



TT 



TT 



--- + 27.-2,1,+ ^^ 



+ 



(3.27) 



_7r^_ll 
4''^ 12 18 



Note that Eq. (3.26) is similar but not identical to the equation derived in Ref. [23] for 
W[Cs] (Eq. 4.5 in Ref. [23]). The equation for W[Cs] takes the form 



d\nW[Cs]{Nfi/m) 



-2A{as{^i)) In^ -r(a,(/x)), 



(3.28) 



where, according to Eq. (3.22), 
Prom Eq. (3.23) we have 



ainTr[C,s](LQ,s(/i)) das 



da. 



din /J, 



\nW[Cs]{l,as{fi)) = CF 



asin) 



TT 



SO the NLO result for r{as{^i)) is: 



r(a,(/i)) = -Cf 



asin) 



TT 



asin) 



TT 



24 



1 vr^ 



Ca f9 



9 ^ 12 ^ /3o U 12 18/ 



(3.29) 



-Cf^ + Cf 



TT 



as{ii) 



TT 



1 



TT 



54 72 



/ 55 



TT 



9 



VIO8 144 4 



3.4 Relation between distribution and fragmentation in the Sudakov region 

In the following we show that the function controlling the large-A^ limit of the pertur- 
bative heavy-quark fragmentation function is identical to S{NiJ,/m) of the heavy-quark 
distribution function. 

According to the discussion in Sec. 3.2 above, the log-enhanced terms in the distri- 
bution function can be computed in the Eikonal approximation, namely using the Wilson 
line definition of Eq. (3.19). Following Sec. 2.3 in Ref. [23] we now show that the same 
applies in the case of the fragmentation function. We then prove that in the Sudakov limit 
the two objects are in fact identical. This is a new result. 

The fragmentation function is defined by [4] (see Eq. (5.2) there): 



d{x;fx) = 



A-2e roo 



2tt 



dy- exp{—ip'^y /x) x 



(3.30) 



^Tr I J2 7+(0| *(0) $*_,(oo, 0) \H{p) + X) {H{p) + X\ %{y, 00) *(y) |0)^ 1 , 
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where dimensional regularization in D = 4 — 2e dimensions is assumed and fi is the renor- 
mahzation scale of the operator. Here the trace is taken over Dirac and color indices, the 
sum is over all hadronic states X that can be produced together with the observed heavy 
hadron H{p), p is the momentum of the latter which is assumed to have no transverse 
component and y is a lightlike vector in the "— " direction. For the Wilson lines we use 
the notation of Eqs. (A.l) and (A. 2), where * stands for complex conjugation; Ref. [4] 
expresses the same Wilson line in terms of transposed color matrices, i.e. 



^'-.(00,0) 



P exp 



-ig I dy A+{y % 

10 



exp (^^9 dy ^tiy )tl 



It is straightforward to check that ^'(O) <t*_y{oo, 0) is gauge invariant if the gauge field does 
not transform at infinity. Since d{x\ jj) has support in x G [0, 1] one can define moments as 
usual, 

/•I 

dxx^~^d{x; /x). 



DN(m;iJ,) = / I 
Jo 



(3.31) 



In the perturbative analogue of d{x; fi), which we denote by dpT{x; ji), H{p) is replaced by 
an on-shell heavy quark h{p) with p^ = m^. Contrary to the non-perturbative definition, 
in perturbation theory h{p) carries color, and so does X. 

Because there are no collinear singularities from the heavy-quark propagator, Sudakov 
logarithms are either due to soft gluons, or to gluons that are collinear with the "— " 
direction. Therefore, as argued in Ref. [23], the Sudakov limit can be safely studied in the 
approximation where the dynamical heavy-quark field in each amplitude is replaced by a 
free quark multiplied by a Wilson line along the quark trajectory: 



*(z) V(^) go) ; *(z) $lp(oo, z) i>{z) , 



(3.32) 



where the path-ordered exponential is defined as in Eq. (A.l). 

Taking the infinite-mass limit implies that the detected on-shell quark h{p) is produced 
from the heavy-quark field in the operator rather than from some vacuum fluctuation. Note 
that when considering the Sudakov limit such fluctuations are irrelevant even if the mass is 
not large, since they necessarily involve gluon splitting into a quark-antiquark pair, which 
is regular at x — > 1. Therefore, in this approximation the free quark field ^ annihilates 
h{p) from the external states. Factoring out the Dirac structure one obtains: 

^Tv|^7+(0|*(0)$l,(oo,0) \h{p) + X){h{p) + X\%{y,^)^{ym)^ 

p+e^P+i/-(0|$;(0,oo)$*_j,(y,0)$l^(oo,y) ^ $;(y, oo) $*_p(oo, y) |0)^ 



p' e 



(0| %{0, oo) ^-y{y, 0) $-p(oo, y) |0)^ 



(3.33) 



where we relied on completeness of the set of states \X) (which close the color trace) and 
on the properties of the Wilson lines in Eq. (A. 3). We therefore find that the x — > 1 
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ll 




y (y"<0) 

Figure 2: Minkowski space-time picture of the Wilson line configuration in the second Hne of 
Eq. (3.34), i.e. W*[Cs]{ip ■ y^J^/^m) = W[Cs]{~ip ■ yn/m), representing the perturbative fragmenta- 
tion function in the infinite-mass Hmit (in the rest frame of the produced quark). The two figures 
describe the situation when y~ is positive (l.h.s) or negative (r.h.s), where path ordering on the 
lightlikc segment I2 from to y corresponds to time-ordering and anti-time-ordering, respectively; 
c/. Fig. 1. 



singular terms in the fragmentation function are summarized by: 

p+dy 



dMx;fi)L^i - fj-^^'^"^''^''^^ (O|$,(0,oo)$_,(y,0)$_^(oo,y)|0)t 
= r e-^^'+^-(i--) (0| %{y, 00) c&,(0, y) $_,(oo, 0) |0)^ 

= /I e^«^^-(i--) (0| $_,(y, 00) $,(0, y) $,(00, 0) |0)^ 

J-00 27r 

= r e-'P^y-^'--^ W[Cs]i-tp ■ y^i/m, a,(M)), (3.34) 

where in the third line we defined q = —p and in the fourth we identified the matrix element 

as VF[C5] that was defined in the context of the distribution function, Eq. (3.19) above. 
Finally, in the last line we returned to the original variable p finding that the function 
W^[(75] is evaluated at —ip ■ yn/m. Since dpT{x; fi) is real, sign inversion of the argument 
of M^[Cs'] is equivalent to its complex conjugation (c/. Eq. (3.20)) so we get 

dpT(x; /x)U_i ~ e^P+i/-(i-) W[Cs]iip ■ yiJi/m, a^if^)), (3.35) 
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the same expression as for /pT(a;;/i) in Eq. (3.18). Therefore, we find that as far as the 
X — > 1 terms are concerned, dp-j-ix; fi) is identical to fpT{x; jj,). Having the relation in the 
second Une Eq. (3.34) above, Ref. [23] has defined 

W[CT]{ip-ytJL/m,as{n)) = (0| $p(y, oo) $j,(0, y) #_p(oo, 0) |0)^. (3.36) 

Here we find that W[Ct] is related to ^[Cs] by complex conjugation: 

W[CT]iip ■ ytx/m,as{ti)) = W[CsY{ip ■ yf^/m,as{f^)). (3.37) 

Finally, converting Eq. (3.35) to moment space ip'^y" — > N (see Eq. (9) in Ref. [30]): 

L>^^(m;/x) = H^{m;iJ,) x W[Cs]{Nfx/m,as{lJ,)) + 0{1/N), (3.38) 

so one obtains the same Sudakov exponent as in the distribution function to any logarithmic 
accuracy. The A^-indcpcndcnt terms are summarized by H^r{m;i^). At 0{as) H^^{m; jj.) 
is equal to H^{m;iJ,) of Eq. (3.24), but this may not persist at higher orders. 

3.5 Comments on the calculation of the Sudakov exponent to NNLO 

In this section we explain how the two-loop coefficient d^^ of Eq. (3.4) was obtained. 
Having established the all-order equality between the Sudakov exponents of the heavy 
quark distribution and fragmentation functions as well as the relation between the definition 
based on a dynamical heavy quark with a finite on-shell mass and the Wilson— line definition, 
there are several way to proceed. We follow two: 

• Perform a two-loop calculation using Wilson lines, as done in Ref. [23] by Korchemsky 
and Marchesini. 

• Extract the non-Abelian^ N — >■ oo singular terms from a recent result for the frag- 
mentation function by Melnikov and Mitov [31], which uses dynamical heavy quarks 
with a finite mass. 

Beginning with the latter, the two-loop calculation of Ref. [31] conveniently suites our 
purpose, as it relies on a process-independent definition of the perturbative fragmentation 
function in dimensional regularization [29], where process dependent power corrections in 
the hard scale (the scale at which the heavy quark is produced) are avoided by taking the 
quasi-collinear limit [29,45]. In this limit the gluon transverse momentum and the quark 
mass are taken small while the ratio between them, which depends on the quark longitudinal 
momentum fraction, is fixed. Ref. [29] established this definition and applied it at 0{as), 
conforming previous result [27] which was obtained from heavy-quark production cross- 
section in e+e" annihilation. Ref. [29] also presented results for Sudakov resummation to 
NLL accuracy. Ref. [30] extended the process-independent calculation of the fragmentation 
function to all orders in the large-/3o limit, which, in particular, fixes the Abelian part of 
df^, see Eq. (3.11) and Table 1 above. 

'^The CpNf term was known already in Ref. [30], and was confirmed by Ref. [31]. 
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Ref. [31] gives a general two-loop result for (ipT(a:;/i) in momentum fraction space. It 
also summarizes in moment space the terms which are non- vanishing at large N in Eq. (65) . 
We note^ that the coupling in this paper is renormalized assuming Nj + 1 dynamical 
massless quarks, where the additional flavor corresponds to the heavy quark. This gives 
rise to CpTjia^ terms which are not accompanied by Nf. Converting to our definition of 
the coupling, with Nf light quarks, these contributions drop out. The remaining terms in 
Eq. (65) of Ref. [31] match the general expression of Eq. (3.6) with Eqs. (3.7) and (3.9) 
above provided that d^^^ is given by Eq. (3.4). 

Here we performed^ a two-loop calculation of VF[C5] of Eq. (3.19), along the lines of 
Ref. [23]. We recall that VF[C5] is a path-ordered exponential with two antiparallel rays in 
the timelike directions ri/^ and — n^u, which are connected by a finite lightlike segment ?/_. 
The timelike direction is determines by the heavy-quark momentum: = mn^, but 
in contrast with the original definitions, Eqs. (2.1) and (3.30), and with the calculation of 
Ref. [31] discussed above, the heavy quark is no more a dynamic field. It is replaced by 
a Wilson line that represents interaction with soft gluons only (collinear singularities do 
arise though from the lightlike segment). As shown in the previous sections VF[C5] captures 
the log-enhanced terms in the heavy quark distribution and fragmentation functions to all 
orders, i.e. to any logarithmic accuracy. Working in D = 4 — 2e dimensions in configuration 
space with the Fcynman rules as in Appendix A of Ref. [23] the standard MS scale fx is 
introduced by making the following replacements: 

n-y — >n-yii, 

^.-^(S'(i-^1' + o(«?)) 

Our final result for the renormalized lniy[Cs] to 0{a1) is summarized by Eq. (3.23) 
above. Given that the description of the calculation in Ref. [23] is clear and detailed, wc 
shall not repeat it here. Instead, we give a full account of the calculation of one diagram, 
diagram 11, the one by which we differ from Ref. [23]. This is done in Appendix D. 

4. Conclusions 

In this paper wc studied pcrturbative aspects of the heavy-quark distribution function, 
which has a central role in precision phenomenology of inclusive i?-decay spectra, primar- 
ily in B — > Xg^ and semileptonic decays. Our results include all-order resummation of 
running-coupling effects, in Eq. (2.9), as well as determination of the Sudakov exponent 
in Eq. (3.2) to two-loop order. The result for d^^ is Eq. (3.4) is now established in two 
entirely different calculation procedures: the one performed here following Ref. [23] using a 
Wilson-line operator in configuration space and the one of Ref. [31] for the fragmentation 

®I wish to thank Matteo Cacciari and Kiril Melnikov for related discussions. 

^As mentioned following Eq. (3.23), the result quoted in Ref. [23] differs from my own in the non-Abelian 
coefficient of the single log term, the one that determines d2- Historically, it was the discrepancy with d2 I 
extracted from Ref. [31] which convinced me to repeat the two-loop calculation of Ref. [23]. I wish to thank 
Gregory Korchemsky for his encouragement. 
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function, using the quasi-collinear limit [29,45] in momentum space. With the jet-function 
anomalous dimension already being known to this order (see e.g. [26]), the Sudakov expo- 
nent in inclusive decay spectra is now determined to the NNLL accuracy. As usual, this 
should be matched by the computation of AT-independent terms at 0{a^), which are not 
yet available. 

In addition to presenting the result for Sudakov resummation in the conventional way 
(Eq. (3.2)) that suites fixcd-logarithmic-accuracy calculations, following previous work on 
DGE [14-16,30] we formulated the resummation in Eq. (3.7) as a scheme-invariant Borel 
sum, where power-like separation between perturbative and non-perturbative contributions 
can be implemented by taking a Principal Value prescription. This approach is particu- 
larly advantageous for inclusive B-decay spectra where non-perturbativc corrections are 
substantial. The application of this approach to phenomenology is already under way [17]. 

The quark distribution function is defined here assuming a finite on-shell quark mass, 
while much of the literature on inclusive decays is based on defining it in the m — > oo 
limit. The evolution properties of these objects are different. Starting in Sec. 3.2 with 
the QCD Sudakov resummation formula for the case of a finite on-shell mass we derived a 
strictly factorized form that is consistent with the m — > oo limit. In the former the Su- 
dakov factor is naturally defined with factorization-scale independent normalization. This 
is realized in Eq. (3.2) and in Eq. (3.7). This is not the case in the strictly-factorizcd 
formula of Eq. (3.15), where the first moment strongly depends on the scale. Strict factor- 
ization (Eq. (3.14)) implies that both the Sudakov factor S and the hard factor H acquire 
double logarithmic dependence on the scale to any order in perturbation theory (see e.g. 
Eq. (3.24)). Therefore, in this case both need to be resummed. 

The most interesting finding of our present investigation is the similarity of the dis- 
tribution and fragmentation functions. First we found that in the large-/3o limit these 
functions are identical. Then, we showed that in the Sudakov limit they are represented 
by the same Wilson-line operator — see Eq. (3.34) — so the Sudakov exponent in the two 
cases is identical to all orders. We emphasize that the diagrammatic realization of this 
relation is non trivial: upon calculating separately virtual corrections to the fragmentation 
process one encounters additional Coulomb-phase contributions in individual diagrams, 
that are absent in the distribution case. Our result implies that these contributions cancel 
out in the sum of all diagrams so they make no effect on the Sudakov exponent. 

In spite of these strong relations the distribution and fragmentation functions, defined 
in a process-independent way in dimensional regularization, are not equal. Their DGLAP 
evolution away from the large-x limit starts differing already at two-loop order. Additional 
differences between these functions appear when considering heavy quark-antiquark pairs 
that were neglected here; these are important in the case of fragmentation for low masses 
and away from the large-x limit. In spite of the similarity in the renormalon structure it 
is hard to imagine that there is any relation between power corrections in the two cases. 
Recall that the distribution function is defined with a single hadron in the initial state and 
a completely inclusive final state, making it a forward hadronic matrix element. On the 
other hand in the case of fragmentation both the detected hadron and the jet are in the final 
state and they interact by exchanging soft gluons throughout the hadronization process. 
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A. Properties of Wilson lines 

We define the Wilson-line operator by 

%izi,Z2) = P exp (^ig dz^A^'iz)^ . (A.l) 

with A^{z) = Aa{z)ta where ta are SU(A^c) generators (tl = ta) in the fundamental rep- 
resentation and Pexp indicates that matrices and fields are path-ordcrcd. The notation 
^p{zi,Z2) assumes^'^ that the direction p is from zi to Z2- Anti-path-ordering is denoted 
by: 

^-p{z2,zi) = Pexp (^-igj^ ' dz^Ai'iz)^ . (A.2) 

The following properties of ^p{zi,Z2) are useful: 

causality ^p{z2, Z3)^p{zi, Z2) = ^p{zi,zs) where Z2 is between zi and z^ 

hermiticity ^l{zi, Z2) = ^-p{z2, zi) (A.3) 
unitarity ^p{zi, Z2)^p{zi, Z2) = 1 . 

B. One- loop integrals with a Borel— modified propagator 

Let us compute the quark distribution on an on-shell heavy quark with a single dressed 
gluon based the definition of Eq. (2.1). First we note that the numerator is identical to the 
one computed in Sec. 3.1 of Ref. [30] with the replacement of A^oA^o by 1^, where n is in 
the "— " direction. Upon using a Borel-modified gluon propagator as well as dimensional 
regularization we get: 

/pT(a:;/^)loK) = ^e""'^'" *(y)^.(0,2/)7-h*(0) \h{p))^ (B.l) 

^ i-oo 47r ' J (27r)^' 

^ + 2 + 2f; (2m^ + fc^)(l + ^) 



(-A;2)i+«((p + k)^ - m?) (-P)i+"((p + kf - w?f 



where the ^ part in the first term and the entire second term originate in the Feynman 



p 



gauge part of the gluon propagator and the rest is specific to the Axial gauge (2.2). 

^°There is some redundancy in the notation as Z2 — zi is parallel to p, but it is, nevertheless, convenient. 
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To perform the momentum integration we use the fact that under the integral, a 
factor of in the numerator becomes (1 — x) while the inverse factor becomes 1/(1 — x 
We therefore have to deal with just one type of integral: 



(B.2) 

Using Feynman paxametrization and shifting the integration momentum q = k + p{l — a) 
and changing the order of integration we obtain: 

I{m\x-aAD)=P^ ^'VI}"t^^ tdaa^-\l-ar' T ^ e^^— )-^-- 



r(6)r(a) 7o ' ' J-oo 27r 

iy- 



(27r)^ {q^ -m?{l-afY+^ 

i(-l)^r(a + 6-f) ^^^^^ _ _ „)6-l(^2(i _ «)2)f -a-6 

(47r)-r(6)r(a) Jo 



(47r)^r(6)r(a) 



(m^) 2 - (B.3) 



where we first performed the momentum integration (observing that the exponential e*^ ^ 
can be replaced by 1 since the result is a scalar). Then we performed the y~ integration 
getting a Dirac 5{a — x). This made the integration over the Feynman parameter a trivial. 
We comment that had we used the Eikonal approximation for the massive propagator, i.e. 
{p + kY — mP' — > 2pk we would have obtained the same answer (same F functions) but the 
dependence on x would have been modified such that x""^ — > 1, not affecting the large-x 
limit. 

Eq. (B.l) can now be computed by the appropriate assignments in Eq. (B.3). The 
result is given by Eq. (2.3). 

C. The splitting function and the cusp anomalous dimension 

The non-singlet splitting function has been recently computed to three loops by Moch, 
Vermaseren and Vogt [33]. In addition, its large-/3o limit is known to all orders and it is 
given by [34] 



„ oo 

7(iV,a) = ^^7n(iV)a"+i=^ 



+ 



^(iV + a) -^(1 + a) (C.l) 
N-1 fa^ + 2a-l 1 {1 + af 1 



1 + a N + a 2 + a N + l + a 



2a 



+ 0(1/(3^ 



where A is the large-/3o limit of the cusp anomalous dimension [23,35,36,38,40] discussed 
below, and a = a{p) = (3cias{/-J-)/TT is the large-/9o coupling in MS. 

The perturbative expansion of the cusp anomalous dimension (the anomalous dimen- 
sion of an operator made of two Wilson lines with a cusp), which is also the large- limit 
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of the quark-quark splitting function, is given by 



Po Jo 



du T(u) 



c^r/^oaTOf) , MS 



Ba{u) 



TT 



TT 



+ 



,(C.2) 



where the coefficient af^ is known from the recent calculation of the splitting function [33] . 
Explicitly, 



3 + f3o 



1 1 



TT 

3 12 



/253 57r2 7^ , ^ 



+- 



192 4 ' ^ ^ ^ 



TT 



11 11 TT^ 

18 18 T'"^'^ 720 



(C.3) 



Ca' 



In the large-/?o hmit [34,40]: 



A{asii^)) = — 



Cf sinna r(4 + 2a) ^n^-l/a2^ 



Po TT 6r(2 + a)2 



+ 0(1/(5^). 



(C.4) 



The Borel representation of A{as{fi)) in the full theory can be written in an expanded 
form as in Eq. (3.12), where c„ represent the terms that are subleading in (3q. Upon 
comparing the latter with the second line in Eq. (C.2) and using Eq. (2.5) we get: 



C2 = ar 



„ „MS _|_ I rMS X„MS 

C3 — 0-3 + 2 + <^2 - > 



(C.5) 



where the term involving S2 = P^^/Pq is due to converting from MS to the 't Hooft scheme 
(see Eq. (27) in [26]). Explicitly this gives: 
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(C.6) 



The presence of a 0{\IPl) term is, of course, a special feature of our specific Borel repre- 
sentation (or of the 't Hooft scheme) . 
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D. Calculation of diagram 11 in Ref. [23] 



Let us focus here on the calculation of the non-Abelian two-loop diagram shown in Fig. 3 
in the Feynman gauge. The calculation is done in configuration space using dimensional 

regularization in D space-time dimensions. The Feynman rules are given in Appendix A in 
Ref. [23]. This diagram involves a triple gluon vertex, where two gluons attach at different 
points zi and Z2 along the time like Wilson line representing the incoming heavy quark 
(with momentum = mn^; = 1) and the third gluon attaches to the lightlike line 
going along the "— " direction at the point z^. 




Figure 3: One of the non-Abelian two-loop diagrams contributing to lnVF[Cs] of Eq. (3.19). It 
corresponds to diagram 11 in Fig. 6 of Ref. [23] (published version). 



We parametrize the points along the Wilson lines as in Ref. [23], z\ = Tin, Z2 = T2n 
and Z2 = T^y, and denote the triple gluon vertex by Z4, obtaining: 

Wn = i^'C.C, J\n [ ir, (.-(„. .)„) ■ (^ - J(^,) (D.l) 



'n 

with 



J{zi) = / d^Z4 n D(zi - Z4) (D.2) 

i=l 
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where the Feynman gauge propagator T){z) is: 
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TiD/2-1) 
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The calculation of J'(zi) using Feynman parametrization is straightforward and yields: 
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with 



A = t(l - t)zf + s(l - s)z| + (s + - s - 

— 2st2;i • Z2 — 2s(l — s — t)z2 ■ z^ — 2f(l — s — t)zi ■ z^ 
= T^t{l - i) - 2stTiT2 + r|s(l - s) - 2n • y(l - s - t)T3 [sr2 + in] . 

Taking the derivative in Eq. (D.l) and scaling ti^2 by — 2n • y one obtains: 
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By first changing the integration variable T2 into cr, where T2 = ctti, and then the integration 
variable ri into to, where 



Tl = UJ 



{I - s-t)T3{as + t) 
a'^s{l - s) - 2sta + t{l - t)' 



one finds that the T3 integration is trivial, giving a factor of 1/(8 — 2D), while the u- 
integration yields: 



/•OO 

/ du;u;3-^(l + a;)2-° 
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The result can be expressed as: 
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where In matches the notation of Ref [23]. We have: In = Kn x Jn with 
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Proceeding with the evaluation of Jn we change the integration variable t into z, 
where t = {I — s)z/{l + z), and then, after changing the order of integration, a into x 
where x = {1 + z){l — a) and s into y where y = s/z, getting: 

/•oo rl+z r^/z , , 6— 2D / \ £)— 4 

Jn=y dzil+zf-^j dxj dyy^/2-2(y(l+2z)-l)(l+y(l-x)j (i+y^l-xf^j . 

(D.9) 

Since the dependence of the integrand on z is simple we perform this integration first. The 
price is having several terms as the z-integration extends between max{0,x — 1} and 1/y. 
Fortunately, most of the terms cancel out by symmetry and we obtain: 

-2 /•°° ri^+y)/y ^ ^ „/ n6-2Z)/ ^\D-4: 

Jii = dy dxy^/2-i^D-3 ^ _ ^ _ ^^2 j p -LQ) 

+ ciy + (l + y(l - x)) (l + y(l - x)^) . 

By changing the integration variable x into z where z = \-\-{V — z^jy and then y into it; 
where y = (1 — z)/^;, the two integrals can be computed exactly, yielding: 
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where the first line is the result of the first integral in Eq. (D.IO) while the remaining terms, 
containing hypergeometric functions, correspond to the second. 

Combining Eqs. (D.ll) and (D.8) we obtain the final result for in, which can be 
readily expanded in e, where = 4 — 2e. The result is: 

ill = — e-^ + (-- + -7E^ e-3 + (-- + ^7r2 + -^l - e'^ (D.12) 

192 V 96 96 ' y V 48 1152 96 '-^ 48 ^ ) ^ ' 

+ V^^^ + 576^ + 144^^ - 24 - 48^^ " 24^^ " 576^ J ^" + ^^"^^^ 
Finally, in terms of the renormalized coupling, Eq. (D.6) takes the form: 

= \caCf iin ■ y/xe^^)^^ In 

= \caCf (^^) ' e^^^ e-2^-^/n, (D.13) 

where in the first line we introduced the MS factorization and renormalization scale accord- 
ing to Eq. (3.39) and in the second we defined L = lnip - yjl/m = InNjl/m with jl = e^^/i, 
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as in Eq. (3.23), knowing that the e-expansion of e ^'^^'^/n is free of 7^; terms. Performing 
this expansion, and subtracting the e — > singular terms, we finally get: 

which differs from Eq. (3.6) in Ref. [23] by the single-log term only. 
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